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6 ■ ABSTRACT 

We report on a limited search for pure-rotational molecular hydrogen emis- 
sion associated with young, pre-main-sequence stars. We looked for H 2 v = 
J = 3^1andJ = 4^2 emission in the mid-infrared using the Texas Echelon- 
Cross-Echelle Spectrograph (TEXES) at NASA's 3m Infrared Telescope Facility. 
The high spectral and spatial resolution of our observations lead to more strin- 
gent limits on narrow line emission close to the source than previously achieved. 
One star, AB Aur, shows a possible (2a) H 2 detection, but further observations 
are required to make a confident statement. Our non-detections suggest that a 
significant fraction, perhaps all, of previously reported H 2 emission towards these 
objects could be extended on scales of 5" or more. 

Subject headings: circumstellar matter — infrared: stars — planetary systems: 
protoplanetary disks — stars: pre-main sequence 
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1. Introduction 

Formation of a circumstellar disk is recognized as a natural step in the process of star 
formation and a vital step toward forming planets. Currently our knowledge of protoplan- 
etary disk properties at anthropically interesting distances of 1-30 AU is relatively limited. 
The best spectroscopic constraints on disk temperature and density structure are dominated 
by material either within 0.5 AU of the central source (Najita et al. 1996; Carr, Mathieu, & 
Najita 2001) or at radii >50 AU (Dutrey, Guilloteau, & Simon 1994). This radial sampling 
results from the tracers utilized to date: dust emission and scattering and CO ro- vibrational 
and rotational emission. Using H 2 rotational lines as a tracer may allow study of disks at 
radii of 1-50 AU. 

Recent work by Thi et al. (1999, 2001a, 2001b) report H 2 v = J = 2 -> and 
J = 3 — > 1 emission from pre-main-sequence and main-sequence stars with circumstellar 
disks. From Infrared Space Observatory (ISO) data, these authors derive gas temperatures 
of 100-200 K and masses up to 2xl0 -3 M . In the case of disks around main-sequence 
stars (Thi et al. 2001a), previously thought to be mostly gas-free (Zuckerman, Forveille, & 
Kastner 1995), the presence of a substantial reservoir of H2 alters ideas on the formation of 
giant planets (Lissauer 2001). Unfortunately, the large aperture and low spectral resolution 
of the observations leave open the question of whether the emission really comes from a disk. 

Three mid-IR H2 pure rotational lines are readily available from high, dry sites if ob- 
served at moderately high spectral resolution: v = J = 3 — > 1 [i.e. S(l)] at 17 ^m, 
J = 4 -> 2 [i.e. S(2)] at 12 /im, and J = 6 -> 4 [i.e. S(4)] at 8 fim. The high spectral 
and spatial resolution possible with ground-based spectroscopy, together with large telescope 
apertures, can result in greater sensitivity to certain gas distributions than satellite obser- 
vations. In particular, ground-based spectroscopy is well suited for detecting point sources 
with narrow line emission, such as gas in Keplerian orbit at a radius of 1-50 AU from a 
solar-mass star. 

We report here on observations of a small sample of young stars taken with the Texas 
Echelon Cross Echelle Spectrograph (TEXES) in an effort to confirm the ISO detections and 
explore the feasibility of ground-based observations of H 2 from circumstellar disks. 

2. Observations 

We used TEXES (Lacy et al. 2002) at the NASA 3m Infrared Telescope Facility (IRTF) 
to observe the H 2 J = 3 — > 1 transition [A = 17.0348 ^m or v = 587.032 cm" 1 ] and the 
J = 4 — > 2 transition [A = 12.2786 /im or v = 814.425 cm -1 ]. All the observations were made 
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with TEXES in high-resolution mode. Pertinent details regarding the observations may be 
found in Table 1. All sources were observed while using the IRTF offset guide camera as 
well as guiding on the dispersed continuum seen through the spectrograph. 

The case of GG Tau deserves special comment since it has weak continuum, is a strong 
case for ISO detection (Thi et al. 1999; Thi et al. 2001b), and has a unique geometry. 
The source is a quadruple system composed of a pair of binaries. GG Tau A is a 0."25 
binary with millimeter continuum emission (Guilloteau, Dutrey, & Simon 1999) and HST 
scattered light observations (Silber et al. 2000) showing a circumbinary ring extending 
roughly from 180 to 260 AU (l."29 - l."86). The ring is tipped at 37°, resulting in an ellipse 
on the sky with the major axis running essentially E-W. When we observed GG Tau A, we 
widened the slit to 3" and rotated the instrument to orient the slit so that the major axis 
of the projected ellipse would lie along our slit. The weak continuum of GG Tau meant we 
could not guide on continuum signal through the spectrograph. Therefore, we repeatedly 
checked our infrared boresight by observations of a Tau, a nearby, bright infrared source. 
The maximum boresight offset found based on the a Tau observations was 0."5, with a mean 
offset of 0."2. By summing over ^4" along the slit, we can confidently state that we observed 
essentially all of the GG Tau A gap region. 

Correction for atmospheric transmission was done using bright infrared continuum ob- 
jects, either stars or asteroids. Although most stars later than spectral type G show pho- 
tospheric features at R> 20, 000, we can locate features using the Kitt Peak Sunspot Atlas 
(Wallace, Livingston, & Bernath 1994) and the ATMOS3 photospheric atlas (Geller 1992) 
and know that there are no features near the H2 transitions. Asteroids have no features at 
our resolution. Both the J = 3 — > 1 and J = 4 — > 2 lines are near telluric atmospheric lines, 
but the Doppler shift of the source, the Earth's motion, and the high spectral resolution 
available with TEXES helped to minimize atmospheric effects. 

Flux calibration was done using a blackbody and standard stars. 

3. Results 

The data were reduced using the TEXES data pipeline (Lacy et al. 2002). We extracted 
spectra from the data with several strategies: optimal extraction of the point source to look 
for emission coincident with the continuum, sums over the nodded region to look for diffuse 
emission covering half the slit length, and selected sums along the slit to look for isolated 
emission offset from the point source. We saw no evidence for extended emission in any 
spectrum, although uniform emission on > 5" scales would not be recovered. 
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Figure 1 presents J = 3 — >■ 1 spectra for the three sources we observed that also have 
reported ISO detections: GG Tau, AB Aur, and HD 163296. To indicate the line flux 
reported by Thi et al. (2001b), we overplot Gaussians with an integrated line flux equal to 
the ISO measurements. Since ISO was unable to resolve the line profiles, we have simply 
assumed Gaussian FWHM matching our resolution (5 or 7.5 km s _1 , depending on slit width) 
and 30 km s _1 , with the Gaussians centered at the systemic Vlsr- The figures clearly show 
a discrepancy between the reported line flux and our observations. 

Of all our spectra, only the AB Aur 12 /xm spectrum shows a possible detection (Fig- 
ure 2). A Gaussian fit to the AB Aur data finds a feature centered at the systemic velocity 
with a FWHM=10.5 km s" 1 and a line flux of 2.0 (±1.0) x 10~ 14 ergs s" 1 cm~ 2 . There 
are no terrestrial lines and no photospheric absorption lines in the comparison star at this 
frequency that might artificially produce an emission feature. We discuss the implications 
of these data, particularly in light of the non-detection of the J = 3 — ► 1 line, below. 

Table 2 summarizes the results for all our observations. In all cases, we measure contin- 
uum levels comparable to past measurements, although the noise for GG Tau is such that 
many pixels must be combined to obtain a significant continuum measurement. We estab- 
lished line flux errors by summing over the number of pixels corresponding to the assumed 
line widths at all parts of the spectrum with comparable atmospheric transmission. We then 
use the line flux upper limit for the more conservative case of a 30 km s _1 FWHM line to 
calculate a maximum mass of H 2 assuming temperatures of 150 K and 300 K. Note that 
these masses assume no extinction and do not include gas in equilibrium with optically thick 
dust. 



4. Discussion 

We observed six sources with reasonable evidence for disks and report no convincing 
detection of H 2 pure-rotational emission. Our non-detections have bearing on the interpre- 
tation of previous ISO reports and on the distribution of gas and dust in these sources. In 
all cases where we set limits, those limits are subject to uncertainties in the assumed gas 
temperature and the assumed line widths. To provide an example of our sensitivity to these 
assumptions, we include two line widths and two gas temperatures in Table 2. 

Three reported ISO detections of H 2 from circumstellar disks, out of three tested, are not 
seen with TEXES. For ISO to detect emission that we would miss, the emission can either 
be extended or it can be spectrally broad. Given that the ISO-SWS aperture at 17 yum 
covered 14" by 27", compared with our effective aperture of 2" by 2" (280 AU by 280 AU at 
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Taurus: d=140 pa), and that ISO detected H 2 rotational lines from molecular clouds with 
even fairly modest incident UV fields (Li et al. 2002), the presence of extended emission is 
possible. Thi et al. (2001b) address this issue using 12 CO 3-2 emission and conclude that 
GG Tau shows no contaminating molecular material, that the bulk of molecular material 
for HD 163296 lies in the disk, but that AB Aur may have a significant contribution from 
an extended envelope. For the case of broad lines, we believe we would have seen lines as 
broad as f^50 km s _1 (see Figure 1 and Table 2). Even for the case of GG Tau, we would 
have seen a line as broad as 30 km s _1 at the 2a level, given the best fit line flux in Thi et al. 
(2001b), although the uncertainties are large enough that these observations may not be in 
conflict, especially for a broad line. Recent H 2 t) = 1 -> J = 3 ^ 1 observations at 2 /jm 
show line widths <15 km s" 1 for GG Tau (J. Bary, in preparation), while CO ro- vibrational 
emission toward AB Aur is unresolved at 12 km s _1 (G. Blake, in preparation). Although 
the relationship between these tracers and the pure rotational H 2 emission depends in detail 
upon the excitation mechanism, these observations show that our assumption of 30 km s -1 
line widths when calculating mass upper limits is likely conservative. As a fiducial point, 
the projected Keplerian velocity at a distance of 1 AU around GG Tau Aa (M = 0.8 M ; 
White et al. 1999) and assuming any circumstellar disk has the same orientation as the 
more extended circumbinary ring (inclination i = 37°), would be v=16 km s -1 . For AB 
Aur, taking M=2.4 M (van den Ancker et al. 1998) and inclination i < 45° (Grady et al. 
1999), the projected Keplerian velocity would be v< 32 km s^ 1 . Of course, emission lines 
broadened by Keplerian rotation will generally have widths broader than their characteristic 
Keplerian velocity. 

Although TEXES is a relatively new instrument, we feel reasonably confident in our 
results. The successful detection of continuum emission from each source at levels consistent 
with past observations means we were pointed at the source and that our flux determinations 
are not dramatically in error. We have detected H 2 rotational emission in sources such as 
Uranus (L. Trafton, private communication) and NGC 7027 (H. Dinerstein, in preparation). 
The ISO detections are the result of special data processing and are near the sensitivity 
limit (Thi et al. 2001b). Because the ISO reports of molecular gas reservoirs around main- 
sequence stars such as (3 Pic (Thi et al. 2001a) have such important ramifications for our 
understanding of planetary formation, the current non-detections in pre-main-sequence disk 
sources strongly argue for follow-up observations of the main-sequence stars observed by ISO, 
all of which are reported to be "medium confidence" detections. Follow-up observations are 
particularly needed for the well studied source f3 Pic, where FUSE observations detect no 
H 2 absorption through the edge-on disk against a stellar O VI emission line (Lecavelier des 
Etangs et al. 2001) and where Na I emission shows a gas disk that seems to coexist with the 
well known dust disk (Olofsson, Liseau, & Brandeker 2001; Heap et al. 2000) 
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In the disk sources we observed, the absence of H 2 emission certainly does not mean 
there is no warm H 2 . Where the dust in the disk is optically thick, an emission line would 
result only if there is spatial and/or temperature separation between the dust and H 2 , such 
as a gap in the dust disk or a hot layer above the optically thick disk (D'Alessio et al. 
1998; Chiang & Goldreich 1997) Gaps are believed to be present in GG Tau and GW Ori 
(Guilloteau et al. 1999; Mathieu, Adams, & Latham 1991); for these sources, the mass upper 
limits constrain the gas within the gap. In general, our observations argue, albeit weakly, for 
little spatial separation or temperature differentiation between the bulk of the H 2 gas and 
dust. 

As described in Section 3, the 12 ^m spectrum of AB Aur has a 2a bump at the 
correct velocity for the J = 4 — > 2 line and with a reasonable width. When taken with 
the 3a upper-limit on J = 3 — > 1 emission from this source, 3.3 x 1CT 14 ergs s" 1 cm -2 , 
the J = 4 -> 2 line flux of 2.0 ±1.0 x 10~ 14 ergs s^ 1 cm- 2 (-la, + la) implies T gas > 
380 K (275 K, 500 K), assuming no differential extinction. This is warmer than the 140 K 
seen in CO ro- vibrational lines (G. Blake, in preparation). For simple thermal emission, we 
would expect the H 2 pure rotational lines to come from cooler regions than CO ro- vibrational 
emission, although 140 K gas would not collisionally excite ro-vibrational CO emission. The 
3a upper-limit on J = 5 — > 3 line flux as observed with ISO, 0.4 x 10~ 14 ergs s _1 cm -2 , 
(Thi et al. 2001b) combined with the possible J = 4 — > 2 detection (—la, +la) suggest 
T<190 K (220 K, 170 K). Given the errors and the resulting inconsistency, we consider 
the AB Aur spectrum to give an upper limit to the J = 4 — > 2 integrated line flux of 
4.0 x 10~ 14 ergs s" 1 cm~ 2 . As well as reobserving the J = 3 — > 1 and J = 4 — > 2 lines, 
ground-based observations of the J = 6 — > 4 line at 8 fj,m, which should be quite strong if 
the gas is truly at 380 K, would help determine if the J = 4 — > 2 is actually as strong as it 
appears. 

While the next infrared satellites will have exceptional mid-infrared sensitivity, they 
may not be able to detect low line-to-continuum ratio features such as may be present in the 
AB Aur 12 /iin data. If AB Aur were observed at R=600, the maximum resolution available 
with the SIRTF InfraRed Spectrograph (Roellig et al. 1998), a SNR of RilOOO would give a 
3a detection of our derived line flux. The recommended mid-infrared imager/spectrograph 
for NGST, with R=1500 (Mather k Stockman 2000), would require SNR of 450 for a 3a 
detection. Higher spectral resolution on NGST, still a design goal, would dramatically ease 
this type of observation. Future ground-based observations with 8-10m class telescopes and 
high resolution spectrographs such as TEXES will improve on the sensitivity described here 
by an order of magnitude (a factor of 100 in time). 

In conclusion, we have attempted to confirm several reports of H 2 pure-rotational emis- 
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sion from young circumstellar disks without success. Our observations generally set more 
stringent limits on gas within ~50 AU of the source and with widths <50 km s _1 than 
obtained with ISO. We have found possible (2a) emission from one source, AB Aur, al- 
though inconsistencies in derived gas temperatures lead us to consider the result as an upper 
limit. The absence of emission suggests the gas and dust have little segregation in tempera- 
ture and/or spatial extent, although more observations are necessary. We have shown that 
ground-based, high-resolution spectroscopy is an effective way to search for narrow, pure 
rotational H 2 line emission such as that expected from circumstellar disks. 

We are happy to acknowledge the day and night staff at the IRTF; without their hard 
work, none of this observing would be possible. We appreciate the comments of Ewine 
van Dishoeck and an anonymous referee. This work was supported by grant TARP 00365- 
0473-1999 from the Texas Advanced Research Program. MJR acknowledges USRA grant 
USRA 8500-98-008 through the SOFIA program. This research has made use of NASA's 
Astrophysics Data System and the SIMBAD database operated at CDS, Strasbourg, France. 
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Table 1. Observation Information 
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a The slit was aligned East- West. Normally, the slit orientation is North- 
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Table 2. Results Summary 
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a From Thi ct al. 2001b with la errors given in parenthesis. 

b Errors given are la per pixel. Continuum level determined from >450 pixels where the atmospheric 
transmission is reasonably good. 

c Upper limits are 3a assuming Gaussian FWHM of 5 km/s and 30 km/s, except as noted. 

d Mass upper limits assume 30 km/s FWHM Gaussian line profile, except as noted, and assume a 
temperature of 150 K or 300 K. 

e Assuming FWHM=7.5 km/s due to 3" slit. 

'Derived from fit to feature at systemic velocity plus 2a (see text). Fitted FWHM=10.5 km/s. 
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Fig. 1. — H 2 J = 3 — >• 1 (A = 17.035 jim) observations of the three sources observed by 
TEXES with reported ISO detections: (a) GG Tau; (b) AB Aur; and (c) HD 163296. In all 
cases, the data (histograms) and overplotted with Gaussians that have integrated line flux 
equal to that reported in Thi et al. (2001b). The narrow Gaussian (dotted) matches our 
instrumental resolution, FWHM=7.5 km/s for GG Tau and FWHM=5 km/s for AB Aur 
and HD 163296. The wide Gaussian (dashed) has FWHM=30 km/s. The Gaussians are 
centered on the systemic velocity. At this wavelength, the spectral order is larger than the 
detector; regions of the spectrum set to Jy are not sampled. The increased noise toward 
negative velocities is due to increasing telluric opacity. 

Fig. 2. — H2 J = 4 — > 2 (A = 12.279 /im) spectrum of AB Aur. One complete order (out 
of 8) is shown, but the region near -40 km/s where the terrestrial atmospheric transmission 
is less than 75% is set to zero. The data (histogram) are overplotted with a Gaussian 
fit. The resulting fit has FWHM=10.5 km s _1 and centroid at the systemic velocity. The 
dotted line shows the spectrum of Capella (divided by 10), the atmospheric calibrator. The 
features in Capella near -80 km/s and +70 km/s are photospheric OH absorption. Although 
we attempted to correct for the OH absorption before dividing AB Aur by Capella, some 
contamination may still be present. 
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